//
// Test Suite for Orientation::isCCW() function
// Ported from JTS junit/algorithm/IsCCWTest.java

// tut
#include <tut/tut.hpp>
// geos
#include <geos/algorithm/Orientation.h>
#include <geos/geom/Polygon.h>
#include <geos/geom/Geometry.h>
#include <geos/geom/CoordinateSequence.h>
#include <geos/geom/Coordinate.h>
#include <geos/io/WKTReader.h>
#include <geos/io/WKBReader.h>
// std
#include <string>
#include <memory>
#include <cassert>
#include <sstream>

using namespace geos::algorithm;

namespace tut {
//
// Test Group
//

struct test_isccw_data {
    typedef std::unique_ptr<geos::geom::Geometry> GeometryPtr;

    std::unique_ptr<geos::geom::CoordinateSequence> cs_;
    geos::io::WKTReader reader_;
    geos::io::WKBReader breader_;

    test_isccw_data()
        : cs_(nullptr)
    {
        assert(nullptr == cs_);
    }

    ~test_isccw_data()
    {
    }
};

typedef test_group<test_isccw_data> group;
typedef group::object object;

group test_isccw_group("geos::algorithm::CGAlgorithms::isCCW");

//
// Test Cases
//

// 1 - Test if coordinates of polygon are counter-clockwise oriented
template<>
template<>
void object::test<1>
()
{
    const std::string wkt("POLYGON ((60 180, 140 240, 140 240, 140 240, 200 180, 120 120, 60 180))");
    GeometryPtr geom(reader_.read(wkt));

    cs_ = geom->getCoordinates();
    bool isCCW = Orientation::isCCW(cs_.get());

    ensure_equals(false, isCCW);
}

// 2 - Test if coordinates of polygon are counter-clockwise oriented
template<>
template<>
void object::test<2>
()
{
    const std::string wkt("POLYGON ((60 180, 140 120, 100 180, 140 240, 60 180))");
    GeometryPtr geom(reader_.read(wkt));

    cs_ = geom->getCoordinates();
    bool isCCW = Orientation::isCCW(cs_.get());

    ensure_equals(true, isCCW);
}

// 3 - Test the same polygon as in test No 2 but with duplicated top point
template<>
template<>
void object::test<3>
()
{
    const std::string wkt("POLYGON ((60 180, 140 120, 100 180, 140 240, 140 240, 60 180))");
    GeometryPtr geom(reader_.read(wkt));

    cs_ = geom->getCoordinates();
    bool isCCW = Orientation::isCCW(cs_.get());

    ensure_equals(true, isCCW);
}

// 4 - Test orientation the narrow (almost collapsed) ring
//     resulting in GEOS during execution of the union described
//     in http://trac.osgeo.org/geos/ticket/398
template<>
template<>
void object::test<4>
()
{
    std::istringstream
    wkt("0102000000040000000000000000000000841D588465963540F56BFB214F0341408F26B714B2971B40F66BFB214F0341408C26B714B2971B400000000000000000841D588465963540");
    GeometryPtr geom(breader_.readHEX(wkt));
    cs_ = geom->getCoordinates();
    bool isCCW = Orientation::isCCW(cs_.get());
    ensure_equals(isCCW, true);
}

// 5 - Test orientation the narrow (almost collapsed) ring
//     resulting in JTS during execution of the union described
//     in http://trac.osgeo.org/geos/ticket/398
template<>
template<>
void object::test<5>
()
{
    std::istringstream
    wkt("0102000000040000000000000000000000841D588465963540F56BFB214F0341408F26B714B2971B40F66BFB214F0341408E26B714B2971B400000000000000000841D588465963540");
    GeometryPtr geom(breader_.readHEX(wkt));
    cs_ = geom->getCoordinates();
    bool isCCW = Orientation::isCCW(cs_.get());
    ensure_equals(isCCW, true);
}

} // namespace tut

